A Modified Block Newton Iteration Forapproximating an Invariant Subspace of Asymmetric Matrixralf L
نویسنده
چکیده
In this paper we propose a Modiied Block Newton Method for approximating an invariant subspace S and the corresponding eigenvalues of a symmetric matrix A. The method generates a sequence of matrices Z (k) which span subspaces S k approximating S. The matrices Z (k) are calculated via a Newton step applied to a special formulation of the block eigenvalue problem for the matrix A, followed by a Rayleigh-Ritz step which also yields the corresponding eigenvalue approximations. We show that for suuciently good initial approximations the subspaces S k converge to S in the sense that sin ' k with ' k := ](S k ; S) converges to zero Q-quadratically provided that the eigenvalues belonging to S are separated from the remaining ones.
منابع مشابه
Certification Using Newton-Invariant Subspaces
For a square system of analytic equations, a Newton-invariant subspace is a set which contains the resulting point of a Newton iteration applied to each point in the subspace. For example, if the equations have real coefficients, then the set of real points form a Newtoninvariant subspace. Starting with any point for which Newton’s method quadratically converges to a solution, this article uses...
متن کاملA New Inexact Inverse Subspace Iteration for Generalized Eigenvalue Problems
In this paper, we represent an inexact inverse subspace iteration method for computing a few eigenpairs of the generalized eigenvalue problem Ax = Bx [Q. Ye and P. Zhang, Inexact inverse subspace iteration for generalized eigenvalue problems, Linear Algebra and its Application, 434 (2011) 1697-1715 ]. In particular, the linear convergence property of the inverse subspace iteration is preserved.
متن کاملA parallel Newton-Krylov flow solver for the Euler equations on multi-block grids
We present a parallel Newton-Krylov algorithm for solving the three-dimensional Euler equations on multi-block structured meshes. The Euler equations are discretized on each block independently using second-order accurate summation-by-parts operators and scalar numerical dissipation. Boundary conditions are imposed and block interfaces are coupled using simultaneous approximation terms (SATs). ...
متن کاملBlock Arnoldi-based methods for large scale discrete-time algebraic Riccati equations
Algebraic Riccati equations (discrete-time or continuous-time) play a fundamental role in many problems in control theory. They arise in linear-quadratic regulator problems, H ∞ or H 2-control, model reduction problems and many others. In this talk we propose numerical methods for large discrete-time algebraic Riccati equations (DARE). We present block projection methods that allow us to comput...
متن کاملA Three-Dimensional Multi-Block Newton-Krylov Flow Solver for the Euler Equations
A three-dimensional multi-block Newton-Krylov flow solver for the Euler equations has been developed for steady aerodynamic flows. The solution is computed through a Jacobian-free inexact-Newton method with an approximate-Newton method for startup. The linear system at each outer iteration is solved using a Generalized Minimal Residual (GMRES) Krylov subspace algorithm. An incomplete lower/uppe...
متن کامل